Cooper pair delocalization in disordered media 
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We discuss the effect of disorder on the coherent propaga- 
tion of the bound state of two attracting particles. It is shown 
that a result analogous to the Anderson theorem for dirty su- 
perconductors is also valid for the Cooper problem, namely, 
that the pair wave function is extended beyond the single- 
particle localization length if the latter is large. A physical 
justification is given in terms of the Thouless block-scaling 
picture of localization. These arguments are supplemented 
by numerical simulations. With increasing disorder we find 
a transition from a regime in which the interaction delocal- 
izes the pair to a regime in which the interaction enhances 
localization. 



I. INTRODUCTION 

It has long been known that the wave function of an 
electron moving in a random potential becomes spatially 
localized. This effect was first predicted by Anderson |^] 
and is termed Anderson localization. In one and two di- 
mensions all quantum states are localized in the presence 
of any amount of disorder while in three dimensions local- 
ization occurs only above some critical disorder. While 
this phenomenon is now fairly well understood in a single- 
particle picture, the inclusion of interactions in the dis- 
ordered many-body system is a non-trivial problem. In- 
sight can be gained by studying the simpler case of just 
two interacting particles in a random potential. In this 
context, it has recently been claimed that the interaction 
can actually lead to a delocalization effect, in the sense 
that the spatial extent of the two-body wave function 
is larger than the single-particle localization length 
This delocalization was found for both attractive and re- 
pulsive interactions, at least for some of the eigenstates 
in the continuous spectrum of energy eigenvalues. Some 
authors objected to this finding j|] and this point is still 
being debated. Furthermore, the relevance of these re- 
sults for the many-body problem has not yet been clari- 
fied. 

In this paper our main concern is the case of an at- 
tractive interaction. The concept of a two-particle bound 
state (the Cooper pair) plays an important role in the 
theory of superconductivity. There is theoretical and ex- 
perimental evidence for the existence of superconductor- 
insulator transitions, where localized states combine co- 
herently into a superconducting condensate with a finite 
superfluid density. This finding motivates the question 
addressed in this paper: will the bound state of two at- 



tracting particles be extended over distances larger than 
the single-particle localization length? We find that a 
result analogous to the Anderson theorem for localized 
superconductors is valid for this problem in the limit 
of large single-particle localization lengths. By increas- 
ing the disorder we find a transition from a regime in 
which the interaction increases the localization length to 
a regime in which the interaction reduces the localization 
length. 



II. LOCALIZED SUPERCONDUCTORS 

The Hamiltonian for a disordered conventional (s- 
wave) superconductor is 



Ho~V V| (r)^| (r) (i-)V'i {r)dr (1) 



H 



where Hq represents the single-electron part including a 
spatially random external potential and the interaction 
parameter V is taken as positive. The attractive inter- 
action could be due to exchange of phonons or purely 
electronic mechanisms but its origin does not concern us 
here. Depending on the disorder, the eigenstates </)n(r) of 
Hq can be extended or localized (with localization length 
Lc at the Fermi level). 

At T=0 the system described by (|l|) is a superconduc- 
tor with spatially constant order parameter A(r) = 
< V'T (i")V'X (r) >= A when the condition 



ANfL;^ > 1 



(2) 



holds m . A is then given by the same expression as that 
for a clean superconductor. This is the Anderson theo- 
rem. The important message of equation (||) is that if the 
disorder is strong enough to localize the single-particle 
states, the superconducting order parameter A and crit- 
ical temperature Tc will remain unaffected as long as 
the number of single-particle states in the energy range 
A contained in the localization volume L"^ is still large. 
(The superfluid density, on the other hand, is greatly re- 
duced.) If the amount of disorder is further increased so 
that (||) is no longer valid the superconducting order pa- 
rameter will fluctuate strongly in space and the critical 
temperature will be lowered. 

Note that because A(r) cx< V'T(i")'0i(r) > is the wave 
fimction of the condensate, the Anderson theorem tells us 
that the attractive interaction is delocalizing the Cooper 
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pairs. Keep in mind, however, that in this derivation 
the Cooper pairs are strongly interacting (and actually 
overlapping) with each other, forming a correlated liquid. 
In BCS theory the concept of Cooper pair only has the 
formal significance that strong pair correlations exist be- 
tween the particles in phase space. In the regime of va- 
lidity of the Anderson theorem A is given by 
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(3) 



The denominator of this expression is the quasi-particle 
excitation energy. The binding energy of an "isolated" 
Cooper pair is given by an expression of the same form 
as this one (see eq. (||) below) but with a different denom- 
inator. This is because equation (^ takes into account 
the interactions between Cooper pairs. 

Experimentally, it has been found that superconduc- 
tivity persists up to the Anderson metal-insulator transi- 
tion ||^-||]. These systems exhibit activated conductivity 
above Tc. The coexistence of superconductivity and lo- 
calization has also been observed in underdoped high-Tc 
superconductors [p^p2[. 



III. 



COOPER PROBLEM IN A RANDOM 
POTENTIAL 



In what follows the problem of two interacting elec- 
trons in a random potential will be addressed for the spe- 
cific case of an attractive interaction. The aim is to see 
how the interplay of disorder and interaction affects the 
coherent propagation of the electrons in the ground state. 
The attractive interaction will be assumed to be short- 
ranged. We can describe this system by an Anderson- 
Hubbard Hamiltonian 



H = -t 



E - 

<r,r'> ,cT 



(4) 



with Cr representing the site energies randomly dis- 
tributed over a width W, and U > 0. The single-particle 
eigenf unctions 0„ are assumed to be localized by the dis- 
order with a localization length Lc. 

We search for the two-particle ground state 



H^iri,r2) = i;o1'(ri,r2) 



(5) 



where ri and r2 are the coordinates of the electrons. Be- 
cause the bound state is a spin singlet, \l/(ri,r2) is a 
symmetric function of its arguments and has the general 
form 



^'(ri,r2) = ^gn,m(l}n{ri)(t)*^ir2) 



(6) 



where (j)n are normalized single-particle eigenfunctions of 
(^ with energy eigenvalues In order to simulate the 
presence of a Fermi sea, it is assumed that the two elec- 
trons can only be paired in states 0„ which lie above the 
Fermi surface. According to the values of the parameters 
in (^) we recognize several different regimes, li U ^ t 
(strong attraction) then the electrons are tightly bound 
and move together like a heavy boson with hopping am- 
plitude t'^/U in an environment with disorder W. This 
boson would then become easily localized. In the remain- 
der of the paper we concentrate on the regime in which 
the interaction is not strong ([/ < t). In this case is 
essentially the result of electron pairing in time-reversed 
single-particle eigenstates of (||) for not too strong disor- 
der. This can be seen as follows. 

The Schrodinger equation for the Cooper pair wave 
function is 

^(ri,r2) 

-En,,n ^ Eo-i„^i,„ 0n(ri)0™(r2). (7) 

It admits the solution 'i'{r, r) — const, if the condition 



Eq — 2^„ 

7V(C,r) 



(8) 



holds. (Here N{^, r) denotes the local density of states 
at the point r). Thus a result analogous to the Anderson 
theorem is valid for the function ^(r,r). In such a case 
the solution to (^ is 



^'(ri,r2) cx^ 



0»(ri)C(i-2) 

Eo — 2^„ 



(9) 



and has no overlap with the state (/)„(ri)(/)^(r2) if n 7^ 
TO. The condition (^) can be satisfied if the integral of the 
local density of states / d^N{£^, r) over an energy interval 
of order | Eq \ does not depend on r. This is possible even 
if the disorder is strong enough to localize the single- 
particle states as long as the condition 



I Eo I NpLj! » 1 

holds. The binding energy | Eq \ will then be the same 
as that for a clean system. So we reach the conclusion 
that the attractive interaction can delocalize the pair or, 
at least, increase its localization length. 



IV. DELOCALIZATION BY INTERACTIONS 

The delocalization by interactions due to the attractive 
interaction can be understood within the block-scaling 
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picture of localization introduced by Thouless ||l3|] . In a 
recent paper , Imry has used the block-scaling picture 
in order to argue that interactions (attractive or repul- 
sive) should, in some cases, delocalize some of the eigen- 
states of the continuous spectrum of the pair of electrons. 
In what follows we extend the argument to the case of a 
Cooper pair. 

Suppose the system is divided into blocks of linear size 
Lc (measured in units of the lattice spacing) so that the 
mean level spacing in a block is Ai = B/ L'j: with B equal 
to the band width. We can then solve the Cooper prob- 
lem, as above, for each block by pairing the two electrons 
in time-reversed states localized inside the block. If the 
fluctuation of the binding energy is smaller than the ef- 
fective scattering amplitude of the pair between blocks 
then the pair will be extended over many blocks. We 
denote the Cooper pair in the j-th block by ^'j. Next, 
we estimate the scattering amplitude tef / , due to the in- 
teraction^ of the Cooper pair between two neighbouring 
blocks as 

W/ = - < *j I uY^-^r^n,,^ I ^-.v > 

r 

= -U^gnngn'n'Y^ I </'n(r) H (r) ^ 

n,n' r 

~ -Yd^dnngn'n'- (10) 

n,n' 

We now note that (;„„ depends smoothly on n and has no 
nodes because it is the wave function of a ground state, 
namely 
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The normalization of implies gnn ~ XL^o | /Ai)^ 
and the number of terms in the sum in (|lO| ) is large, of 
the order of (| Eq \ /Ai)^. All those terms interfere 
constructively yielding 

\Eo\U_ U\E^\ 
" Ai Ld ~ B • 

The fluctuation of the binding energy | Eq \ due to the 
randomness in ^„ is of the order of 6E ~ y/\ Eq \ Ai (see 
the Appendix). Thus the condition for delocalization is 
obtained as 

teff U 

The ratio U/B is the BCS product NpV. The following 
points should be noted: 

(i) There are no phase correlations to consider in the 
sum ( p^ because all the terms are real and negative; 

(ii) The large effective hopping amplitude ie// resulted 
from gnn being a smooth function of n with no nodes. In 



other words, because we have been considering the ground 
state. 

(Hi) If the sign of the interaction U is reversed, to make 
it repulsive, we do not expect any delocalization to oc- 
cur in the ground state. The reason for this is that '^j 
for each block (using time-reversed pairing) would now 
essentially involve only the single-particle eigenstate 0o 
with the lowest energy: ^'j(ri,r2) « 0o(ri)0o(''2)- Then 
teff ^ C^/Lc ~ Ai if [/ « S. If we remove the constraint 
of time-reversed pairing then Imry's argument would pre- 
dict delocalization only if the pair has a certain excitation 
energy above the Fermi level ||l^,|l^ . So it would not be 
in the ground state. 

V. NUMERICAL RESULTS 

We have also performed numerical calculations of the 
Cooper pair wave function for ID systems of up to i = 
100 sites using the Lanczos algorithm. In order to im- 
pose the constraint of a filled Fermi sea the single-particle 
eigenstates for a given realization of the disorder poten- 
tial are calculated and the eigenenergies of the states be- 
low the Fermi energy are shifted by a large amount. Thus 
only the one-particle states above the Fermi surface are 
accessible for the Cooper pair. Before doing the Lanczos 
diagonalization of the two-particle problem the Hamil- 
tonian is transformed back to real space representation 
where the number of nonzero matrix elements is much 
smaller than in the basis of the eigenstates of the non- 
interacting disordered system. 

In order to determine the spatial extent of the Cooper 
pair wave function with respect to both the relative and 
the center-of-mass coordinates we have calculated two 
different quantities. The first one is the participation 
ratio 

1 

'■"EJ*(r,r) |4 

where ^'(r, r) is normalized to unity. Since only the diag- 
onal part (ri = r2) of the wave function is involved in the 
calculation the participation ratio can also be interpreted 
as the localization length of the Cooper pair. 

The second quantity, which is related to the relative 
coordinate of the two electrons, is the average size d of 
the Cooper pair defined as 

d=J2 I*(ri,r2) nri-r^l. 

ri,r2 

Our results are obtained for a system of 100 sites and 
each data point represents the average over 200 realiza- 
tions of the disorder potential. Since the fluctuations of 
the participation ratio r are very large we found it more 
convenient to average log r instead of r itself. 

Fig. 1 shows the behavior of logr as a function of U/t 
for different values of the disorder strength W/t. While 
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FIG. 1. Participation ratio r of the diagonal component of 
the Cooper pair wave function as a function oiU /t for several 
values of the disorder strength W/t. 
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FIG. 2. Average size d of the Cooper pair as a function of 
U/t for W/t = 1, 1.5, 2, 2.5 (from top to bottom). 



in the case of strong disorder the interaction leads to a 
decrease of the participation ratio, for smaller values of 
W/t an enhancement of logr is observed, at least for not 
too large U/t. This is in qualitative agreement with our 
analytical arguments. One should however keep in mind 
that the calculations are done for a rather small system 
and that finite size efi^ects can be of importance, espe- 
cially in the weak disorder region where the localization 
length becomes large. 

In Fig. 2 we see that the size of the Cooper pair de- 
creases rapidly with increasing interaction strength. The 
values for U/t ^ which should be of the order of the 
one-particle localization length are considerably reduced 
due to finite size effects. 



VI. CONCLUSION 

The interplay of disorder and interaction is a complex 
problem and therefore it is useful to start with the two- 
body problem in a disordered medium. The effect of a re- 



pulsive interaction seems to be very different from that of 
an attractive interaction since the latter can induce prop- 
agation of the bound state of two particles (their ground 
state). The binding energy and spatial extent of the 
bound state is insensitive to disorder beyond the point 
where single-particle states become localized. This result 
for the Cooper pair corresponds to the familiar Anderson 
theorem for the many-body problem of dirty supercon- 
ductors and is valid in a regime where the attraction 
and the disorder are not too strong so that electrons are 
paired in time-reversed single-particle eigenstates. These 
localized states can then be combined coherently into an 
extended pair wave function (eq.(^). In analogy with 
the case of dirty superconductors discussed in y] where 
the superfiuid density is greatly reduced by the disorder, 
we also expect the energy of the wave function (|) to be 
much less sensitive to the boundary conditions than that 
in a clean system. In other words, disorder has a stronger 
effect on the sensitivity to boundary conditions than on 
the localization length. It is possible to go away from 
this "Anderson regime" in two ways. If the interaction 
is increased then both the binding energy and the effec- 
tive mass of the pair increase, so the pair becomes more 
localized. On the other hand, if the interaction is kept 
not too strong but the disorder is increased, the binding 
energy and the localization length are reduced. 
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Appendix 



We want to prove that SE ~ ^/EAi. 
energy E is obtained from the equation 
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U 



E 
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E + 2e„ 



E + 2^' 



The binding 



(11) 



The density of states is N{^) = No{0 + n{0 where iVo(0 
is the average density of states and n(^) the noise due to 
disorder. Then E = Eq + SE where Eq is the binding 
energy derived from iVo(^). Expanding the right-hand 
side of the above equation for small 6E and taking into 
account the definition of Eq we arrive at 



SE 



n{0 
Eo + 2^ 



(12) 



We make the following assumptions about the moments 
(averages) of n(^): < n(^) >= and < n(^)ri(^') >= 
noexp(— (^ — ^')^/A^). This latter condition is only in- 
tended to express the fact that the correlation persists 
over an energy of the order of the single-particle level 
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spacing. Taking the square of (H) and averaging we ob- 
tain 



This gives SE - f^y/A^. The quantity no/Np should 
not be larger than unity. 
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